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Abstract 

We describe a search method for fast moving (/3 > 5 x 10~^) magnetic monopoles 
using simultaneously the scintillator, streamer tube and track-etch subdetectors 
of the MACRO apparatus. The first two subdetectors are used primarily for the 
identification of candidates while the track-etch one is used as the final tool for 
their rejection or confirmation. Using this technique, a first sample of more than 
two years of data has been analyzed without any evidence of a magnetic monopole. 
We set a 90% CL upper limit to the local monopole flux of 1.5 x 10~^^ cm~'^ s^^ sr~^ 
in the velocity range 5 x 10~^ < /3 < 0.99 and for nucleon decay catalysis cross 
section smaller than ~ 1 mh. 
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1 Introduction 



Within the framework of Grand Unified Theories (GUT), supermassive mag- 
netic monopoles (m ~ lO^^GeV) would have been produced in the early 
Universe as intrinsically stable topological defects when the symmetry of the 
unified fundamental interactions was spontaneously broken [1]. At our epoch 
they should be searched for in the cosmic radiation as remnants of primordial 
phase transition (s). The velocity range in which GUT monopoles should be 
sought spreads over several decades [2]. If sufficiently heavy (m > lO^^GeV), 
GUT monopoles will be gravitationally bound to the galaxy with a velocity 
distribution peaked at /9 ~ 10^^. Lighter monopoles (m < 10^^ GeV) would 
be accelerated in one or more regions of coherent galactic magnetic field up 
to velocities oi f3 > 10"^ while other acceleration mechanisms (e.g. in the 
neighborhood of a neutron star) could bring them to relativistic velocities. 

MACRO was a multipurpose underground detector (located in the Hall B of 
the Laboratori Nazionali del Gran Sasso, Italy) optimized for the search for 
GUT monopoles with velocity P = v/c > 4 ■ 10~^ and with a sensitivity below 
the Parker bound (i.e. 10~^^ cm~^s~^sr~^ [1)2]). The apparatus was arranged 
in a modular structure with overall dimensions of 76.5 x 12 x 9.3 m^ and was 
made up by three subdetectors: liquid scintillation counters, limited streamer 
tubes and nuclear track detectors (CR39 and Lexan) [3]. 

In this work we describe an analysis technique which uses all of the MACRO 
subdetectors simultaneously in the search for fast moving magnetic monopoles 
(i.e. /3 > 5 X 10~^) in the cosmic radiation. The results on a first data sam- 
ple is also reported. The analysis procedure uses the data coming from the 
scintillator and the streamer tube subdetectors to identify candidate events. 
This is done by reconstructing the energy release (using scintillators' Energy 
Reconstruction Processor -ERP- as well as streamer tubes' Charge and Time 
Processor, QTP) and the particle's trajectory (using the streamer tubes' digi- 
tal hit information). The selected tracks are then searched for in the track-etch 
layers as a final tool for their rejection or confirmation. 

Since the three techniques are independent and complementary, we obtain a 
good rejection power against the background due to high energy muons and 
a high reliability of the possible candidates. In the following sections we will 
describe the signatures of a fast monopole in the MACRO detector and the 
analysis procedures. Analysis efficiencies and acceptance evaluations will then 
be given together with the results of a search performed on a first sample of 
more than two years of data. 
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2 Fast monopole signatures in MACRO 

A complete description of the MACRO detector, with particular emphasis on 
magnetic monopole detection techniques, has been given elsewhere [3-5] . Here 
we will concentrate on the aspects which concern this particular analysis. 

The energy loss suffered by a fast magnetic monopole in matter is due to 
atomic excitation and ionization. Thus the average rate of energy loss fol- 
lows immediately from an extension of the Bethe-Bloch formula for a moving 
magnetic charge. This allows a straightforward and well-grounded calcula- 
tion of the detector response. The peculiar characteristic of the fast magnetic 
monopoles is their large ionizing power compared either to considerably slower 
monopoles or to minimum ionizing electrically charged particles. The energy 
released in a MACRO streamer tube by a /3 > 5 • 10"^ monopole is at least one 
hundred times larger than that due to a minimum ionizing particle (m.i.p.) 
[6]. Since the charge of the streamer pulse has been shown to have a logarith- 
mic dependence on the ionizing power [6,7], it provides a good rejection of the 
muon background. The response of the scintillation counters to a magnetic 
monopole of a given velocity is reported in [8,9]. The expected signal for a 
monopole is about a factor 30 greater than the m.i.p. level for /9 = 5 ■ 10~^ and 
is even larger for greater velocities. In the case of the track-etch detectors, the 
calculated value of Restricted Energy Loss (REL, defined as the fraction of 
the particle's energy loss contained in a 10 nm diameter around its trajectory) 
for a magnetic monopole with p — 5 ■ 10~^ is ~ 60 MeV cm^g~^, and it in- 
creases rapidly at higher velocities [10]. This value is well above the detection 
threshold of the MACRO CR39 nuclear track detector, REL„i„ = 26 MeV 
cm^g~^. The Lexan has a much higher threshold compared to that of CR39, 
making it sensitive to relativistic monopoles only. 

The particular signatures of a fast monopole in MACRO are thus: a) a single 
track in the detector's spatial and temporal views corresponding to a particle 
with a velocity greater than 5 • 10~^ c, and b) a large energy deposition 
along the particle trajectory. 

The background for this search is almost entirely due to the fiux of high 
energy single muons. At the depth of MACRO, about 5% of these have en- 
ergies greater than the threshold energy above which the discrete energy loss 
processes of nuclear interactions, bremsstrahlung, and pair production start to 
dominate, viz., ~ 1 TeV. These discrete mechanisms produce showering events 
in the apparatus. 

This background is generally reduced by a combination of geometrical and 
energy cuts imposed on each of the subdetectors. An analysis based on all three 
subsystems in MACRO allows the use of these cuts in a rather conservative 
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way. In particular, the energy cut remains safely within the linear response 
of the system. Moreover, any systematic errors are greatly reduced by the 
combination of measurements from the three subsystems. 

Energetic muons {Ef^ > 1 TeV) showering in the apparatus are efficiently re- 
jected by requiring a clean single track in the streamer tube system. However, 
this cut will also reject monopole events in which electromagnetic or hadronic 
showers are produced. In particular, a monopole with 7 > 10 (i.e. /? > 0.99) 
can produce a 5-recy in a single electron collision with enough energy to produce 
an electromagnetic shower, giving rise to secondary tracks. Furthermore, pro- 
cesses such as hadronic interactions and e"*" — e~ pair production are effective 
for 7 > 100 [11] and will produce secondary tracks. Finally, another possible 
source of non-clean monopole-induced single tracks in the apparatus might 
be due to the nucleon decay catalyzed by monopoles following the Rubakov- 
Callan effect [1] . Due to these effects the present analysis might lose efficiency 
in identifying monopole events if the catalysis cross section acat > 1 or 
p > 0.99. 



3 Event reconstruction 

This search begins with event selection performed by both the ERF and the 
Streamer tube Horizontal Monopole Trigger (SHMT). 

Event reconstruction starts from the procedures used in the search for slow 
monopoles with the streamer system [4]. Once event tracking has been per- 
formed, the energy loss observed in the scintillators and the analog response of 
the streamer tubes are used to achieve the necessary rejection power against 
the muon background. 

3.1 Event tracking in the Streamer Tube System 

The details of the event tracking are fully described in [4]. The strategy 
adopted to eliminate the radioactivity background is briefly summarized in 
the following. The SHMT selects events with at least 7 horizontal planes with 
a rough temporal alignment (a "/? slice"). For each of the triggered /3-slices, 
wire hits are selected on the basis of time measurements made by the QTPs, 
while strip hits are identified by requiring a spatial match with the wire view. 
If a track is found in each of the two spatial views, a more refined time track 
is searched for by using the QTP timing information from the selected spatial 
hits. Furthermore we require a minimum of 7,6 and 7 points along the track 
in the wire, strip and time views respectively. Only events satisfying these re- 
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quirements and consisting of a single track reconstructed in each of the wire, 
strip, and time views are analyzed further. 



3.2 Measurement of the energy loss in scintillation counters 

The ERP, a system measuring the energy deposition in each MACRO single 
scintillator and time of flight across the apparatus, has been extensively de- 
scribed elsewhere [3,5]. In order to obtain AE, the energy lost by a particle in 
a scintillation counter, the raw ADC values from each of the two tank ends are 
corrected for light attenuation in the scintillator and a normalized weighted 
mean of the two is taken. The normalization constant for each scintillator is 
updated on a weekly basis by fitting the Landau-like signal size distribution 
produced by a selected set of cosmic muon events. Correcting for the light 
attenuation requires knowledge of the longitudinal position of an event in a 
scintillator. This can be obtained either from the streamer tube tracking sys- 
tem or from the difference of arrival times of the light at the two scintillator 
ends. In this analysis, we use the former since it yields better energy reso- 
lution than the time-difference method. In Fig.l the difference between the 
positions reconstructed with these two different methods is shown for a sam- 
ple of cosmic muons. The width of this distribution is primarily due to the 
time resolution of the scintillators and it yields at ~ 750 ps. In Fig. 2 the value 
of the reconstructed energy deposition of a sample of muons is shown as a 
function of the longitudinal position of the event in the scintillator. As can 
be seen the correction for varying sensitivity as a function of this position is 
performed correctly in almost all the sensitive volume. Near the scintillator 
ends this procedure fails by underestimating the energy lost by the particle. 
For this reason a fiducial volume is defined by cutting the last 20 cm close to 
each end. 

A study of the response of the ERP to large signals has been fully described 
in [5]. The response of the whole chain from the photomultiplier tubes (PMT) 
to the recording analog-to-digital converting (ADC) electronics was measured 
as a function of light level by pulsing each individual counter with UV light 
coming from a nitrogen laser (A ~ 337 nm) after passing through a precision 
light attenuator. The system has been shown to have a linear response up 
to a level corresponding to roughly 5 times the amount of light produced in 
a pathlength of ~ 30 cm by a minimum ionizing muon. In this search the 
synergetic use of the digital (tracking) and analog (pulse charge) information 
coming from the streamer tube system allowed putting analysis cuts below 
the onset of ERP non-linearity. 

Since in this search we are interested in monopoles with velocity down to 
P — 5 ■ 10~^, the scintillator crossing time might exceed the width of the 
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ERP integration gate. In this case only a fraction of the PMT pulse would 
be integrated, thus underestimating the energy lost in the scintillator. The 
ERP ADC integration gate is 200 ns wide and the time spent by the digital 
electronics to produce it is 70 ns. As a consequence the PMT signal to be 
integrated is delayed by 75 ns. Moreover the decay time of the slow component 
of the scintillation light, the electronic chain bandwidth, and the geometry of 
the detector actually set a minimum pulse duration of about 60 ns, which can 
be estimated by observing the muon pulses. This implies that saturation effects 
might start for pulses produced by particles which traverse the scintillator in 
more than 135 ns. However in 135 ns, a (3 — 5 ■ 10~^ monopole travels about 
20 cm and loses at least ~ 465 MeV [8]. Since the analysis cut will be set 
at a much lower level, the width of the integration gate does not affect the 
detection efficiency. 

The selection criteria adopted in this analysis for the ERP hits are then sum- 
marized as follows: 

• The values given by the ERP TDC's at the two different system thresholds 
{high and low) cannot differ, at each scintillator end, by more than 3 ns. 
This is essentially a check that the system is working properly. 

• The measurements of the longitudinal position of a hit as performed by the 
streamer system tracking and by the ERP TDC's, cannot differ by more 
than 100 cm (see Fig.l). 

• A fiducial volume is defined to guarantee a proper energy reconstruction. 
This is done by excluding 20 cm close to the scintillator ends (see Fig. 2). 

• The pathlength of the particle in the scintillator volume, as measured by 
using the streamer system track parameters, must be greater than 10 cm. 

For all the events in which at least two different scintillators are selected, the 
time of flight information is also available. 

3.3 Measurement of the streamer pulse charge 

3.3.1 Study of the analog response of the streamer tubes 

Experimental studies of the dependence of the streamer charge on the par- 
ticle's trajectory and ionizing power arc reported in [6,7,12]. In this section 
we show the results of a study of the MACRO streamer tubes charge, as 
measured by the QTP system, using cosmic ray muons. This is done to bet- 
ter understand the detector characteristics as well as the dependence of the 
charge on the geometry of the ionizing track, the gas mixture parameters, etc. 
The understanding of these parameters allowed us to set up a procedure for 
distinguishing a m.i.p. from a highly ionizing particles such as a fast GUT 
monopole. 
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For each selected event, QTP hits along the track are identified in order to 
measure the streamer pulse charge. They must: 

• spatially match with the wire and strips hits; 

• be in the time view within ±2 us from the temporal track; 

• not be preceded in the previous 5 fis by another hit in the same QTP channel 
(this is required in order to both reject afterpulses and avoid systematic 
errors in the charge measurements due to the ADC's recovery time). 

In the analysis we considered only wire hits with cluster size CLS=1 or 2 in a 
given QTP channel. In the following we will refer to the streamer pulse charge 
as measured by the QTP receiving the signals from the horizontal tubes, once 
the above mentioned selection criteria are applied. 

For slanted particle trajectories with respect to the direction of the anode 
wire, the probability to produce more than one streamer in the gas volume 
increases. As a consequence, after a region (which we refer to as the dead 
zone) in which the multi-streamer production is strongly reduced by space 
charge effects due to the first avalanche, the pulse charge linearly increases for 
increasing Lp, the length of the projection of the particle track (in the sensitive 
volume) along the direction of the wire. This dependence has been studied in 
reference [12] using a muon beam. Here we report also the results of a study 
on the MACRO streamer tubes performed with cosmic muons crossing the 
detector. In particular, for each 0.5 cm interval in Lp, the charge distributions 
produced by cosmic muons have been fit. Gaussian fits have been separately 
applied for CLS=1 and CLS=2 hits. The mean value of the fits are given as 
a function of Lp in Fig.3. These behaviors are in good agreement with the 
aforementioned results obtained in [12]. The dependence of the charge on Lp 
comes out to be: 

Q^Qo if Lp<Ld (1) 

Q^Qo + S-(Lp- La) if Lp > L^ 

where Qo, S and L^ are the values of the charge produced by vertical muons, 
the slope of the linear part of the plot, and the width of the dead zone, respec- 
tively. Identical behaviors are obtained with vertical streamer tubes. Only the 
results for Lp < 5 cm are reported. This is due to the angular distribution of 
the cosmic muons that limits the statistics for Lp > 5 cm. It must be noticed 
that, given the trigger requirements of the SHMT, the restriction to the anal- 
ysis of events with Lp < 5 cm results only in a 6% reduction of the geometric 
acceptance. 

The procedure used to obtain the plots shown in Fig.3 has been applied to 
the whole data sample. For calibration purposes the values of Qo, S and L^ 
have been measured on a run-by-run basis. This allowed us to study and mon- 
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itor the MACRO streamer tube response, its dependence on the gas mixture 
parameters and the correlation between the CLS=1 and CLS=2 values. The 
behaviors of Qq, S and L^, for CLS=1 and CLS=2, are reported in Fig. 4 
as a function of the solar time. The fluctuations in the value of the charge 
produced by vertical muons are due to small variations of the parameters of 
the gas mixture. In March 1994, during hardware work for detector upgrade, 
there was an abnormal streamer charge response (see Fig.4), due to a strong 
air contamination of the gas mixture. This caused the tubes to operate in 
Geiger mode. 

Small variations in the actual values of Qoj S and are very useful in order 
to study their mutual dependencies and to compare them with the results 
obtained in the literature [12]. One example is the dependence of S on Qq, 
which is shown in Fig. 5 for CLS=1 hits. Since the increase of the pulse charge 
for larger values of Lp is due to the production of more than one streamer 
avalanche along the particle trajectory, one should expect S ~ Qo/Ld- How- 
ever, in agreement with what is reported in [12], this approximation fails for 
large values of Qq, when it is not easy to take into account space charge effects 
in a simple phenomenological model. The important feature of the first plot 
of Fig. 5 is that for a highly ionizing particle (i.e. large Qq) the value of S has 
to be larger than the corresponding value for a m.i.p. in the same detector 
configuration. This is particularly important for this analysis since it allows a 
better discrimination between fast monopoles and cosmic muons. 



3.3.2 Measurement of the average charge along the track 

The results of the study of the MACRO streamer tube response and the afore- 
mentioned calibration procedure allowed us to correct the pulse charge mea- 
surements for both geometrical effects and gas mixture variations. 

The charge measured by the QTP system for the selected hits has been nor- 
malized by dividing it by the value of the average charge produced by a muon 
in the same operational and geometrical conditions. The muon charge is calcu- 
lated by using Eq. 1 with the value of Lp for the event under consideration and 
the values of Qo, S and from the cahbrations for a given CLS and run. As 
can be seen in Fig.4, the values of Qq, S and for CLS=2 have large spreads 
around their mean values. This is a purely statistical effect in the fit to the 
various charge distributions. However the ratios of the CLS=2 constants to 
the CLS=1 ones are very stable in the whole data set. As a consequence, we 
avoided small effects due to the aforementioned spreads by using the values 
obtained for CLS=1 multiplied by the measured scahng factor. 

In order to measure the streamer tube response to the passage of a particle, 
we defined a variable F = Q/Q^, where Q is the charge actually measured for 
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a given event and is that produced by a muon and calculated as described 
above. As was reported in [6,7] F has a logarithmic dependence on the particle 
ionizing power. In particular, for a monopole with /3 > 5 • 10~^, F > 4 is 
expected. The value of F has been measured in each selected event by using the 
QTP information coming from the horizontal streamer tubes for CLS=1 and 
CLS=2 hits Q. In this way n different measurements for CLS=1 (F^^^^^^) and 
/ different measurements for CLS=2 {^cls=2) have been obtained along the 
particle trajectory. The value of F associated with the event is then obtained by 
averaging these n + l measurements, weighted by the inverse of the variances of 
the charge distributions, a^^^g^^ and (JcLS=2y obtained from the data sample 
under consideration. In order to reduce the fluctuations of the value of F, we 
do not include in the average the hit with the largest value of Tcls=sI cycLs=s 
(s=l or 2). Therefore 

Y _ Sr=l ^hLS=ll ^CLS=l + S!=l ^hLS=2l ^CLS=2 ~ CLS=s / <^CLS=s)raax , 

which is calculated on np =n + l — l different points. Given the considerations 
concerning the efficiency and the acceptance of the analysis (see next sections), 
we required an event to have at least rir = 6. 

The distribution of F for a data sample is shown in Fig.6. As expected, the 
shape is almost Gaussian with a tail towards the largest values that is due to 
both fluctuations and 5 ray production by high energy muons. As anticipated, 
a magnetic monopole with (3 >b ■ 10"'^ would have F > 4. However, detailed 
studies of the streamer formation process and of the effects of the read-out 
electronics on the pulse shape show that the actual value of the cut must be 
a suitable function Vcut{Lp) as discussed in the following. 



4 Analysis flow 

The analysis method can be summarized as a sequence of three different 
phases: run selection, event selection and reconstruction, and candidate se- 
lection. 

In the first phase all runs longer than 0.5 hours are selected (the average time 
length of a MACRO run is 6.5 hours). Evident electronic problems are ruled 
out by discarding runs with abnormal trigger rates. This selection results in a 
total live time of more than 80 % of the solar time. 



^ For CLS >2 the poor statistics (due to the muon angular distribution) does not 
allow setting up a calibration procedure on a run by run basis. 
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Event selection requires that both the SHMT and the ERP trigger have fired. 
Then single tracks must be identified in the spatial and temporal views with 
the requirements listed in Sec. 3.1. Given a space track, at least one scintillator 
must satisfy the selection criteria discussed in Sec.3.2. For each of the selected 
ERP hits, the energy reconstruction procedure is applied. If at least two differ- 
ent scintillators fired, a precise time of flight measurement (at the level of less 
than l.ns) is also possible. The QTP hits for the streamer pulse charge mea- 
surements are then selected (see Sec.3.3.1). This allows the measurement of 
^CLS=s CLS=1 and CLS=2 wire hits in the horizontal streamer tubes, 

by using the streamer charge calibration constants previously produced for 
each run. The calculation of the F of the event is then made by making an 
average on at least rir = Q QTP hits (see Sec. 3. 3. 2). The average amount of 
events reconstructed following these requirements is ~ 3 ■ lO^evt/yr. 

Monopole candidate events are then selected. The energy reconstructed in each 
selected scintillator counter must be AE > {AE)min = 150 MeV (see Sec. 5. 3). 
On the ~ lO^evt/yr which are left a cut on F is imposed: F > Fcut(Lp) (see 
Sec. 5. 4). Finally the analysis of the possible survived candidates (~ 3evt/yr) 
is made by the scanning of the track-etch sheets identified by the streamer 
track. 



5 Efficiency evaluations 

As can be envisaged from the previous sections, muon events pass the first two 
analysis steps but not the last, which is based on energy loss measurements 
in the scintillator, streamer and track-etch subdetectors. This allowed the use 
of the cosmic muons as an important tool for an evaluation of the analysis 
efficiency. The data from [7,12] and some characteristics of the streamer tube 
read-out card have been used to compute the efficiency of the cut on the analog 
streamer response. 

5.1 Streamer tube system: trigger and event reconstruction efficiencies 

The SHMT sensitivity to muons allows a simple and direct measure of the 
overall reconstruction efficiency Egt of the analysis steps based on the streamer 
information. It has been estimated according to the procedure fully discussed 
in [4]. Basically Sgt is expressed as the product of the efficiency of the recon- 
struction algorithm Srec times the efficiency Strg due to the trigger and all the 
electronics used to produce the streamer data: Sgt = Srec ■ ^trg- Actually the 
two contributions are not independent, leading to a conservative estimate of 
the overall efficiency. The reconstruction efficiency Srec for particle with /3 ~ 1 
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is measured by the ratio between the rate of muons reconstructed by the anal- 
ysis and the expected value. Then a check is performed of the independence 
of Srec on the velocity of the particle. This is done by allowing the natural ra- 
dioactivity hits to simulate time alignements corresponding to slow particles. 
Electronic failures that can be responsible for trigger inefficiency, i.e. Strg, are 
monitored by computing the ratio between the measured trigger rate and the 
expected one on the basis of the knowledge of the trigger circuitry and of the 
radioactivity background rate. 



5.2 The ERP efficiency 



The scintillator counters considered in a given event are identified by using the 
tracking information provided by the streamer tubes. The single scintillator 
efficiency Sscim with respect to the event reconstruction procedure with the 
streamer tubes, is measured for each run and each counter by checking if the 
ERP trigger associated with the tanks crossed by the reconstructed streamer 
track actually fired. This single scintillator efficiency, Ssdm is very close to 
one except for a few counters with evident hardware problems. The ERP 
efficiency Eerp, has then been evaluated for each run by measuring the fraction 
of scintillators with Esdn > 95%. This procedure is very sensitive to the overall 
ERP efficiency and it automatically takes into account all the inefficiencies of 
the tracking algorithm. 



5.3 The ERP energy cut 



As calculated in [8], a magnetic monopole with /3 > 5 • 10~^ produces a light 
yield in a scintillator which is more than 30 times larger than that due to the 
passage of a m.i.p. If we take into account that this estimate has a factor of 
two uncertainty [8], the minimum light yield produced by a monopole with 
the shortest accepted pathlength in the fiducial scintillator volume (10 cm) is 

MeV 

AE^vf > 15 X 1.8 — — -■ psc ■ 10 cm > 232 MeV (3) 
g/cm^ 



where Psc — 0.86 g/cm^ is the mass density of the scintillating liquid. For 
this value of AE, the non linearity effects (due to the PMT and/or ADC 
saturation) are below the 20% level and the energy resolution is about 5% [5]. 
Therefore a cut which imposes AE > 150 MeV has an efficiency greater than 
99%. 
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5-4 The streamer charge cut 



As shown before, a monopole with ~ 5 • 10~^ would produce a track in 
the streamer system with F ~ 4. The streamer charge would be even larger 
for faster monopoles. The analysis requires F > TcutiLp). The suitable value 
of Fcut(Lp) is chosen to ensure an efficiency of at least 90% for recognizing 
a monopole. To evaluate Tcut{Lp), a complete and detailed simulation has 
been performed of the full analysis chain starting from the production of the 
electronic signal. To ensure that the input waveforms for the simulations were 
correct, we collected 80 MHz WFD data, at different Lp values, from a streamer 
tube operating in a muon beam [12] under conditions very similar to those in 
MACRO. 

Since a) a monopole with /3 ~ 5 • 10~^ is expected to produce streamer pulses 
with a charge four times higher than that produced by a muon, b) the dead 
region length is essentially independent on Qq, and c) the slope S increases 
with Qq, the monopole charge parameters are conservatively set at: 



L^-L>^ (4) 



In order to simulate the waveform of the signal produced by a monopole, the 
muon streamer response was multiplied by the ratio of the charge produced 
by a monopole and that by a muon as computed with Eq.l and Eq.4. In the 
simulation, the effect of the protection diodes of the read-out cards is also 
considered. They actually introduce an amplitude saturation of wire signal 
outputs corresponding to more than three times the signal produced by a 
m.i.p. (see Fig. 7). For each Lp value both muon and monopole waveforms 
have been conservatively saturated at a level corresponding to two times the 
amplitude of the single streamer pulse. 

The values of F have been determined, separately for muon and monopole 
events, by using the charges of the simulated pulses at different Lp. For this 
purpose a detailed Monte Carlo simulation of the MACRO detector provided 
the number of hits with CLS=1 and CLS=2 for muons and monopoles, taking 
into account the different angular distributions and the different effects of the 
trigger requirements for monopoles and muons. 

A result of this procedure for /3 = 5 • 10~^ monopoles is shown in Fig. 8. The F 

distributions of muons and monopoles are well separated for small Lp while the 
effect of the saturation due to diodes is clearly visible for larger Lp. Starting 
from these data, a function Tcut{Lp) was extracted for which at least 90% of 
monopoles have F > Tcut- For monopoles with /3 > 5 ■ 10~^, a larger value 
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Period. 


Runs 


^EHP 


^st 


live \^^^ "-A^ ^ ) 


11 Dec 92 - 13 Jun 93 


5520 - 6329 


0.911 


0.788 


3456 


12 Dec 94 - 24 Jul 95 


8907 - 10556 


0.915 


0.894 


3738 


24 Jul 95 - 23 Sep 96 


10557 - 12816 


0.915 


0.946 


8816 



Table 1 

Main characteristics of the analyzed data sample. The values of the analysis effi- 
ciencies (see Sec. 5) and the integrated live times are also reported. 

of r is expected. However, in the whole j3 range covered by this analysis, we 
conservatively used the Tfyut{Lp) estimated for /3 = 5 ■ 10~^ as shown above. 



6 Analysis results on a first data sample 

The analysis procedure described above has been performed on a first data 
sample collected since December 1992 to September 1996, with the exclusion 
of the period June 1993-December 1994, when the efficiency of the streamer 
and scintillator subdetectors were affected by the works due to the upgrading 
of the apparatus with the construction of the upper part. In particular, by 
considering the ERP and SHMT rates, the ERP efficiency and the calibration 
of the streamer charge response we have selected three data taking periods 
with a total live time of more than 16,000 hours. The main characteristics 
of the considered data set, as well as the values of the analysis efficiencies 
discussed above, are hsted in Tab.l. 

Once all the selected events were fully reconstructed, the 150 MeV energy 
cut mentioned above was applied. Then the streamer charge on the surviving 
events was analyzed. The result is shown in Fig. 9 where the dependence of F 
on Lp is reported together with the function Tcut previuosly described. In the 
first period, three events are close to the cut, and, in the third period, three 
others survive the cut. For these six events the track-etch wagons identified by 
using the streamer tracking system were extracted. The passage of a magnetic 
monopole would cause a structural damage. Chemical etching should result in 
the formation of collinear etch-pit cones of equal size on both faces of each foil, 
or of a through-hole, after a prolonged etching. Fig. 10 shows microphotographs 
of holes obtained by exposing a CR39 foil to relativistic sulphur ions (REL 
~ 300 MeV cm^g~^) impinging normally to the detector surface. A candidate 
track must satisfy a three-fold coincidence of the position, incidence angles 
among the layers and should also give the same value of REL. A detailed 
description of the analysis procedure is given in [13]. For each extracted wagon 
the top layer of CR39 was etched in NaOH 8N at 80° C for ~ 120 hours. No 
track compatible with the crossing of a monopole was found. Since the detector 
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threshold is at more than 9 standard deviations below the expected REL for 
a monopolc with /3 ~ 2 • 10~^, the efficiency of the track-etch analysis can be 
considered 100% in this search. 

As a further check, the measured values of energy loss (in the scintillator 
counters) and [3 of surviving candidates have been compared with the expected 
monopole signal as a function of its velocity [8,9]. This comparison is reported 
in Fig. 11. For five of the six events the Time of Flight (T.o.F.) information 
was provided by the streamer system alone, since only one fired scintillator 
counter was present. In this case the error on the reconstructed velocity is large 
because of the limited time resolution of the streamer tubes. As can be seen for 
all the events the measured energy loss values are well below the expectations 
for monopoles. The six events are interpreted as high energy muons producing 
energetic 5 rays. 



6.1 Flux upper limit calculation 

To evaluate the acceptance, 30 000 events have been simulated and written 
with the format of real data in order to use the same off-line code employed 
for the analysis. Since in the first data set the upper part of MACRO {at- 
tico) was not in operation, configurations both with and without it have been 
simulated. By taking into account all the cuts required by the analysis, we 
find an acceptance of ~ 3565 m^sr. In Fig. 12 the acceptance dependence on 
the minimum imposed rir is reported. In the analysis > 6 is required. By 
considering all the effects due to electronics, triggers, cuts on AE and F, etc., 
the product of the efficiency times the live time integrated in the three periods 
is 11853 hours. Since no candidate was found, the monopole flux upper limit 
at 90% C.L. is 

< 1.5 ■ 10-^^ cm-^s-^r-^ (5) 

for 5 • 10~^ < /9 < 0.99 and monopole induced nucleon decay catalysis cross 
section a < 1 mb. For /3 > 0.99 and a > 1 mb the efficiency of this search 
might be reduced by an amount which is difficult to estimate due to theoretical 
uncertainties (see Sec. 2). 

This is the first search in the literature performed with three different detection 
techniques. Furthermore, in the velocity region 10~^ ^ /3 < 0.99, the limit is 
more stringent than any other obtained with scintillation or gas detectors (see 
for instance [14-17]). 
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7 Conclusions 

A combined analysis technique to search for fast magnetic monopoles with 
the MACRO detector has been described. The use of scintillator, streamer 
tubes, and track-ctch data at the same time, and the large energy release of 
fast monopoles, provide a clear signature that allows us to reject the cosmic 
muon background with highly reliable and safe cuts. A first sample of more 
than two years of data has been analyzed. No candidate survived thus setting 
a local monopole flux upper hmit at 90% CL of < 1.5 • 10~^^ cm~^s~^sr~^ 
for 5 • 10~^ < /5 < 0.99 and nucleon decay catalysis cross section a < Imb. 
Even with this reduced data sample in the region 10^^ ^ /3 < 0.99, the limit 
is more stringent than any other obtained with scintillation or gas detectors. 
By using all of the MACRO data (i.e. up to December 2000), the described 
analysis technique will be able to extend its sensitivity to fluxes of the order 
of 5 ■ 10^^^ cm~'^s~^sr~^. 
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Fig. 1. Difference between the measurements of the position of the particle crossing 
along the scintillator as given by the ERP TDC and by the streamer tube tracking 
information. The time resolution of the scintillators comes out to be at c:^ 750 ps. In 
the analysis the hits for which this difference is greater than 100 cm are not used. 
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Fig. 2. Result of the energy reconstruction procedure in the sensitive volume. The 
reconstructed energy is independent on the position of the particle crossing along the 
scintillator longitudinal axis. On the right, the cut applied to exclude corner clipping 
tracks and border effects is also shown. 
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Fig. 3. Charge produced by cosmic muons in the MACRO horizontal streamer tubes 
as a function of Lp. The two plots refer to CLS=1 and CLS =2 data. 
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Fig. 4. Qo, S and as a function of time, for CLS=1 (black points) and CLS=2 
(grey points) hits. The black strips on the top indicate the periods used for this 
search. 
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Fig. 5. Dependence of S on Qo for CLS=1 hits in the horizontal tubes. This behavior 
is also compared with the ratio Qq/L^. 
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Fig. 6. Distribution of the values of F for a data sample. A monopole with 
(3 > 5 ■ 10~^ would have F > 4. No cuts on the energy deposition in the scintil- 
lators has been applied at this stage. 
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Fig. 7. The effect of the protection diodes of the read-out card is shown. The output 
amplitude saturates at a level corresponding to three times the pulse produced by a 
muon. 
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Fig. 8. r values for simulated muons and /3 = 5 - 10^^ monopoles (in grey) for three 
different Lp values. The effect of the saturation of the pulse amplitude introduced by 
protection diodes (see text) is clearly visible. 
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Fig. 10. Microphotographs of holes in CR39 from relativistic S^^~^ ions after 120 h 
of etching in NaOH 8N at 8(F C. The sulphur holes are clearly distinguishable from 
the surface background. 
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Fig. 11. Light yield of the six selected candidates compared with the expected 
monopole signal as calculated in [8,9]. As can he seen for all the events the measured 
values are well below the expectations for monopoles. 
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Fig. 12. Fraction of the effective geometrical acceptance after all cuts as a function 
of the minimum required nr. Stars refer to data simulated without considering the 
attico. The acceptance for np > 6 is 3565m?' sr. 
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